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ABSTRACT 

An  information-geometric  approach  to  sensor  management  is  in¬ 
troduced  that  is  based  on  following  geodesic  curves  in  a  manifold 
of  possible  sensor  configurations.  This  perspective  arises  by  ob¬ 
serving  that,  given  a  parameter  estimation  problem  to  be  addressed 
through  management  of  sensor  assets,  any  particular  sensor  config¬ 
uration  corresponds  to  a  Riemannian  metric  on  the  parameter  man¬ 
ifold.  With  this  perspective,  managing  sensors  involves  navigation 
on  the  space  of  all  Riemannian  metrics  on  the  parameter  manifold, 
which  is  itself  a  Riemannian  manifold.  Existing  work  assumes  the 
metric  on  the  parameter  manifold  is  one  that,  in  statistical  terms,  cor¬ 
responds  to  a  Jeffreys  prior  on  the  parameter  to  be  estimated.  It  is 
observed  that  informative  priors,  as  arise  in  sensor  management,  can 
also  be  accommodated.  Given  an  initial  sensor  configuration,  the  tra¬ 
jectory  along  which  to  move  in  sensor  configuration  space  to  gather 
most  information  is  seen  to  be  locally  defined  by  the  geodesic  struc¬ 
ture  of  this  manifold.  Further,  divergences  based  on  Fisher  and  Shan¬ 
non  information  lead  to  the  same  Riemannian  metric  and  geodesics. 

Index  Terms —  Information  geometry;  Sensor  management 

1.  INTRODUCTION 

The  work  of  Amari  and  others  [1]  on  the  use  of  methods  of  Rie¬ 
mannian  geometry  to  analyze  statistical  estimation  problems  is  of 
increasing  interest  to  researchers  in  signal  processing.  This  method¬ 
ology,  known  as  information  geometry ,  provides  a  rigorous  frame¬ 
work  for  measuring  the  power  of  data  to  discriminate  values  of  pa¬ 
rameters.  These  ideas  date  back  to  Rao  [9],  who  showed  that  the 
Fisher  information  of  a  likelihood  used  in  an  estimation  problem  can 
be  seen  as  a  Riemannian  metric  on  the  parameter  manifold. 

This  paper  brings  an  information-geometric  perspective  to  a  class 
of  sensor  management  problems  by  casting  the  objective  of  sensor 
management  as  parameter  estimation  and  describing  how  this  leads 
to  the  role  of  sensor  management  as  selecting  a  Riemannian  metric 
for  the  parameter  manifold.  Established  results  in  Riemannian  ge¬ 
ometry  [3],  outside  the  context  of  information  geometry,  show  that 
the  collection  of  all  Riemannian  metrics  on  a  Riemannian  manifold 
is  itself  an  (infinite-dimensional)  Riemannian  manifold.  In  problems 
where  the  collection  of  possible  sensor  actions  is  suitably  modeled 
by  a  smooth  finite-dimensional  manifold,  the  space  of  interest  is 
a  finite-dimensional  sub-manifold  of  this  infinite-dimensional  Rie¬ 
mannian  manifold.  A  perspective  is  developed  in  which  the  best 
sensor  management  action  to  take,  in  terms  of  gathering  the  most 
information  relevant  to  the  estimation  objective,  is  locally  charac¬ 
terised  in  terms  of  geodesic  curves  in  this  space. 

Much  of  the  development  in  subsequent  sections  of  this  paper  is 
rather  abstract  and  draws  upon  mathematical  machinery  that  is  unfa¬ 
miliar  to  many  researchers  in  sensor  management  area.  To  provide  a 


more  concrete  context  in  which  to  illustrate  some  of  the  concepts  that 
arise  in  later  sections,  is  it  helpful  to  begin  by  setting  forth  an  exam¬ 
ple  problem.  Suppose  two  mobile  sensor  platforms  and  one  station¬ 
ary  target  (emitter)  are  located  in  the  plane  R2,  as  depicted  in  Fig.  1. 
The  goal  is  to  estimate  the  position  of  the  target  from  bearings-only 
measurements  taken  at  the  sensors.  Since  the  sensors  are  mobile,  the 
sensor  management  problem  is  to  identify  the  trajectories  of  sensor 
motion  that  will  yield  the  best  estimate  of  the  target  position.  More 
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Fig.  1.  An  illustrative  scenario  involves  estimating  the  position 
(xe,ye)  of  a  single  stationary  emitter  from  bearings-only  measure¬ 
ments  received  at  two  mobile  sensing  platforms  located  at  respective 
coordinates  (xi,  j/i)  and  (x2,  yf)  in  the  plane. 

specifically,  the  target  position  is  ( xe ,  ye )  and  the  sensor  positions 
are  ( Xj,yj )  for  j  =  1,  2.  Denoting  Xj  =  Xj  —  xe  and  yj  =  yj  —  ye, 
the  bearing  of  the  target  from  sensor  j  is  ipj  =  arctan(i),, / xf).  The 
sensor  measurements  are  independent  and  von  Mises  distributed, 
each  with  common  inverse  dispersion  parameter  k  and  with  the  cir¬ 
cular  mean  of  the  measurement  at  sensor  j  having  circular  mean  < pj . 

The  following  sections  proceed  to  describe  the  nature  of  a  sen¬ 
sor  model  from  an  information-geometric  viewpoint,  to  define  the 
parameter  manifold,  the  sensor  manifold,  its  metric  structure,  and 
to  derive  a  differential  equation  that  characterizes  geodesic  curves 
on  the  sensor  manifold.  The  development  departs  from  the  purely 
geometric  treatment  in  [3]  in  that  it  allows  for  informative  prior  dis¬ 
tributions  on  the  parameter  manifold  rather  than  restricting  attention 
to  a  volume  form  corresponding  to  the  Jeffreys  prior.  Further,  the 
Riemannian  metric  with  respect  to  which  geodesics  maximize  “en¬ 
ergy  integrals”  on  the  sensor  manifold  is  shown  to  arise  from  both 
Kullback-Leibler  and  mutual  information  perspectives.  Throughout 
this  process,  the  example  just  introduced  will  be  used  to  illustrate 
these  concepts  in  a  concrete  fashion. 


2.  SENSOR  MODEL 

Consider  the  problem  of  estimating  a  parameter  9  from  data  x  col¬ 
lected  by  sensors  up  to  time  t.  Beginning  with  a  prior  probability 
distribution  for  9,  which  may  reflect  what  is  known  from  previous 
measurements  or  side  information,  the  effect  of  taking  measurement 
at  time  t  is  to  provide  a  posterior  probability  distribution,  which 
will  be  assumed  to  be  represented  by  a  posterior  probability  den¬ 
sity  p(9 \x).  If  the  option  exists  to  use  one  of  a  parametrized  set  of 
sensors  or  sensor  configurations,  each  of  these  will  produce  its  own 
posterior  density.  When  these  posteriors  are  known,  selecting  a  sen¬ 
sor  configuration  to  use  for  a  measurement  amounts  to  choosing  one 
of  them  from  which  to  estimate  9.  The  parameter  9  will  henceforth 
be  assumed  to  be  an  element  of  an  m-dimensional  smooth  manifold 
M  (C°°  will  be  assumed,  although  C2  is  sufficient  for  most  of  the 
discussion  here),  which  will  be  called  the  “parameter  manifold.”  In 
the  example  problem,  M  is  R2  because  the  parameter  6  =  (xe,ye) 
is  a  physical  location  in  the  x  —  y  plane. 

Denoting  by  l  =  log  p(x\9)  the  log-likelihood  for  this  problem, 
the  Fisher  information  is 

•Fe  =  Ep(.|0  )[dol®del\,  (1) 

where  dgl  denotes  the  derivative  of  l  with  respect  to  the  parameter 
9.  This  is  well  known  to  be  equivalently  expressed  as 


Fe  =  -Ep(.|0)[V@f]. 


In  this  expression,  V@  represents  the  covariant  derivative  along  any 
connection  in  M;  all  choices  of  connection  give  the  same  quantity 
as  (1).  F$  is  always  a  non-negative  definite  m  x  m  matrix,  and  in 
what  follows  it  will  be  assumed  to  be  non-singular,  thus  giving  rise 
to  a  continuous  family  of  inner  products  over  the  tangent  spaces  of 
the  manifold.  Direct  calculation  in  the  example  problem  shows  that, 
in  the  coordinate  system  described  above  and  depicted  in  Fig.  1,  Fg 
has  the  form 
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where  R2  =  x2  +  y2. 

Through  this  mechanism,  the  choice  of  a  particular  sensor  leads 
to  the  association  of  a  positive  definite  matrix  with  each  9  £  M, 
thereby  imbuing  M  with  a  Riemannian  metric  that  measures  the 
ability  of  that  sensor’s  data,  at  least  locally,  to  discriminate  between 
parameter  values.  It  is  possible  to  calculate  the  shortest  distance,  in 
terms  of  this  metric,  between  two  values  9.  As  discussed  in  [1],  the 
Kullback-Leibler  divergence  between  p{x\Q)  and  p(x\9')  is  approx¬ 
imately  half  of  the  square  of  the  distance  between  9  and  9'  when  9 
and  9'  are  close. 


in  [3],  A  point  in  AT  is  a  Riemannian  metric  on  M\  i.e.,  it  asso¬ 
ciates  a  positive  definite  form  gg  with  each  9  £  M.  Under  suitable 
assumptions,  a  metric  on  AT  is  defined  by 

Gg(h,k)=[  Tr(g^1hggg1kg)  vol(gg),  (3) 

J  M 

where  vol(pe)  =  det(</6>)  d9.  Specific  assumptions  guaranteeing 

finiteness  of  this  integral  are  beyond  the  scope  of  this  discussion,  and 
it  will  suffice  for  the  purposes  here  to  assume  directly  that  it  is  finite. 

Although  the  nature  of  AT  appears  formidable,  realistic  sensor 
management  problems  do  not  require  one  to  work  with  this  entire 
space,  rather  with  a  finite-dimensional  sub-manifold  that  inherits  the 
metric  (3)  from  AT .  The  assumption  that  leads  to  this  situation  is  that 
the  collection  of  all  possible  sensor  configurations  is  parametrized  by 
a  smooth  manifold  S,  which  will  be  called  the  “sensor  manifold.”  In 
the  example  problem,  the  sensor  configuration  is  completely  speci¬ 
fied  by  the  positions  of  the  two  sensor  platforms  in  the  plane;  i.e., 
by  a  =  (*i,2/i,  *2, 3/2)  £  R4.  In  this  case,  the  sensor  manifold 
is  S  =  R4  and  the  only  elements  of  AT  of  relevance  are  those 
metrics  on  M  that  arise  from  a  sensor  configuration  a  in  this  four¬ 
dimensional  manifold. 

Beginning  with  a  sensor  configuration  a  £  S  gives  rise  first  to 
a  likelihood  pa  (x\ 9)  and  consequently  to  a  Riemannian  metric  g(a) 
on  the  parameter  manifold  M,  as  described  in  Section  2.  As  a  Rie¬ 
mannian  metric  on  M,  g(a)  is  an  element  of  AT  (M).  This  mapping 
g  :  S  — t  AT  taking  er  to  g(a)  will  be  called  the  “sensor  geometry.” 
In  what  follows,  g  will  be  assumed  to  be  smooth  and  one-to-one  and 
g(S)  a  sub-manifold  of  AT.  Weaker  assumptions  are  possible  (e.g., 
g  is  an  immersion),  but  full  generality  is  is  not  needed  here  to  ade¬ 
quately  illustrate  the  method.  Through  the  sensor  geometry  map,  the 
finite-dimensional  manifold  S  inherits  the  Riemannian  structure  of 
AT;  i.e.,  the  distance  between  two  sensor  configurations  in  a\  and 
(j 2  in  S  is  taken  to  be  the  distance  between  g{o\)  and  g(a 2)  in  AT. 
This  construction  endows  the  sensor  manifold  S  with  its  own  Rie¬ 
mannian  metric  which  captures,  in  information-theoretic  terms,  the 
“complementariness”  of  sensor  configurations. 
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Fig.  2.  The  sensor  geometry  map  g  allows  the  finite-dimension  man¬ 
ifold  S  of  sensor  configurations  to  inherit  a  Riemannian  metric  from 
the  infinite-dimensional  manifold  AT  of  all  Riemannian  metrics  on 
the  parameter  manifold  M. 


3.  THE  SENSOR  MANIFOLD 

It  has  been  shown  [3]  that  the  collection  AT  (M)  of  all  Riemannian 
metrics  on  the  manifold  M  is  an  infinite-dimensional  (weak)  Rie¬ 
mannian  manifold.  The  structure  of  its  tangent  space  is  described 


4.  GEODESICS 

The  objective  of  determining  good  trajectories  in  the  sensor  mani¬ 
fold  S  will  be  addressed  by  relating  these  to  geodesic  curves  in  AT. 
Following  [3],  consider  a  smooth  curve  7  :  [0, 1]  — >  AT.  For  each 


t  G  [0, 1],  7 (t)  is  a  Riemannian  metric  on  the  parameter  manifold  M 
and  thus  associates  a  positive  definite  matrix  7(f)g  with  each  point 
in  9  G  M.  The  energy  integral  along  the  curve  7  is 

£7  =  \  [  f  Tr(7-177-17))  dF(6)dt.  (4) 

A  JO  J M 

In  this  expression,  dF(-)  is  a  probability  density  on  M ,  7  means 
7 (t)g,  and  7  is  the  derivative  of  7  with  respect  to  t.  Geodesics  in  M 
minimize  E-y,  and  a  variational  approach  is  used  in  [3]  to  obtain  the 
differential  equation  7  =  77_17  for  7 (t),  which  implies 

7(t)  =  7(0)  exp(7(0)_17(0)f) 

The  right-hand  side  of  this  differential  equation  is  observed  to  be  a 
Christoffel  symbol. 

The  induced  metric  at  a  point  a  G  S  is 

Ga{u,v)=  [  Tr  (g(ay1g*(u)g(a)~1g*(v))  dF{9) 

J  M 

where  u  and  v  are  in  the  tangent  space  TSa  of  S  at  a  and  g *  is  the 
push-forward  of  g  :  S  — >  A4.  For  a  smooth  curve  7  :  [0, 1]  — >  S, 
the  energy  integral  restricts  to 

Ei  =  \  f  [  Tr  (fl(7(t))~1g*(7(t))<?(7(t))-1ff*(7(£)))  dF(9)dt 
Z  JO  J  M 


and 

:F»=X'd^:FeXi+X'Wi  BVi' 


d  j.  _  kA(k)  /  0  sirupj  \ 
dxj  R3  Vsin^'  —2  cos  <pj  J 

d  _  kA(k)  /— 2  sirupj  cosi^A 
dyj  9  ~  R3  V  cos  Vi  0  J 


Fig.  4  shows  trajectories  obtained  for  a  particular  case  of  the 
example  scenario.  The  target  is  stationary  at  (1,1),  and  the  sensors’ 
prior  distribution  on  the  target  location  is  normal  with  mean  (1,1)  and 
covariance  0.011.  Sensor  1  starts  at  (0,1)  and  Sensor  2  starts  at  (1,0), 
and  initial  directions  of  motion  are  defined  by  the  geodesic  for  this 
configuration.  The  sensors  move  in  this  direction  for  a  fixed  period 
of  time,  a  new  set  of  directions  is  determined  from  geodesic  calcu¬ 
lations  based  on  the  new  configuration,  the  sensors  move  again,  and 
so  forth.  The  dotted  trajectories  are  extrapolations;  they  indicate  the 
directions  defined  by  the  geodesic  computation  at  the  last  iteration 
computed. 


The  geodesics,  which  are  the  extremal  curves  of  Ey,  satisfy 

7  =  “A  (7, 7) 

where  r  denotes  the  Christoffel  symbol  for  the  Levi-Civita  connec¬ 
tion  on  S. 

In  terms  of  local  coordinates  in  S,  geodesic  equation  in  S  may  be 
obtained  by  solving  a  variational  problem  on  the  path  u.  To  set  this 
up,  it  is  convenient  to  abuse  notation  and  define  a  smooth  function 
u  :  [0,  l]2  — >  S  such  that  u(s,t)\s=o  =  u(t).  With  this  notation, 
and  coordinatizing  Ga  as  Qij, 
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If  this  expression  is  set  to  zero,  further  algebraic  simplification  leads 
to  a  differential  equation  (in  coordinates)  that  characterizes  geodesics 
in  S : 


Fig.  3.  Sensor  trajectories  based  on  geodesic  approximation  for  the 
example  scenario.  Sensor  1  starts  at  (0,1)  and  Sensor  2  starts  at  (1,0). 
The  target  is  at  (1,1). 


Returning  to  the  example  pictured  in  Fig.  1,  the  local  coordinates 
in  S  =  R4  are  xi,  yi,  22,  and  1/2-  The  positive  definite  matrix  g{u) 
corresponds  to  the  Fisher  information  matrix  Tg  given  in  (2).  The 
inverses  and  derivatives  needed  are  calculable,  with 

jr-1  =  (  RlR2  ^  x 

\k71(k)  sin2(y>r  -  p2) ) 

—  sin2  <pj  sin  ipj  cos  pj  \ 
sin  pj  cos  pj  —  cos2  pj  J 


5.  DIVERGENCES  ON  M 

The  proposed  scheme  for  sensor  management  involves  following,  at 
least  locally,  geodesic  curves  in  S  defined  by  the  Riemannian  metric 
S  inherits  from  Ai.  Geodesics  maximize  energy  integrals  of  the 
form  (4),  so  it  is  desirable  to  understand  how  optimization  in  this 
sense  relates  to  the  amount  of  information  gathered  by  the  sensor. 
Consider  first  the  Kullback-Leibler  divergence  D  (.A/2 1 1 Ax )  for  two 
multivariate  normal  distributions  with  equal  means  and  respective 
non-singular  covariance  matrices  g  and  h.  This  is  given  by 

iTr(st-‘— I)  +  ilogj7 


J 


where  |  ■  |  denotes  determinant  and  I  is  the  identity  matrix.  A  diver¬ 
gence  on  M  may  be  defined  by 


^kl(<?,  h) 


dF{6). 


Here,  the  two  positive  definite  matrices  g  and  h  are  regarded  as  aris¬ 
ing  at  each  point  of  M  from  two  Riemannian  metrics.  It  is  evident 
that  ziKL(p,  g)  =  Q  ?mAdgAKL(g,h)\h=g  =  dhAKL(g,h)\h=g  =  0. 
The  corresponding  Riemannian  metric  on  M  is 


dgAK  lU=9  =  dlAKL\h=g  =  \  f  Tr  (g  1  g  g  1g)  dF(0), 

Z  J  M 

as  appears  in  (4). 

Similarly,  one  can  define  a  divergence  on  M  motivated  by  mu¬ 
tual  information  by 


sensor  management  applications.  Navigation  along  geodesic  curves 
in  a  Riemannian  manifold  maximizes  an  energy  integral  involving 
the  metric.  We  have  constructed  two  distinct  divergences  on  M  cor¬ 
responding  to  familiar  information-theoretic  quantities  (Kullback- 
Leibler  divergence  and  mutual  information)  that  have  been  used  by 
various  authors  as  criteria  in  designing  sensor  scheduling  algorithms. 
Both  of  these  are  shown  to  lead  to  the  same  Riemannian  metric  on 
M,  suggesting  the  information  gathering  merit  of  sensor  schedul¬ 
ing  based  on  following  geodesic  curves  defined  with  respect  to  this 
metric. 

While  the  work  presented  here  is  mostly  conceptual,  we  have 
shown  enough  specifics  of  how  the  proposed  method  manifests  in  a 
concrete  example  to  indicate  its  feasibility.  We  are  continuing  to  de¬ 
velop  complete  application  examples  while  simultaneously  working 
out  rigorous  specifics  of  some  of  the  mathematical  foundations. 
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h ) 


1(1  + g^h) 


+  log 


(I  +  h-'g) 


dF(e ) 


This  “symmetrized”  mutual  information  expression  is  equivalent  to 


AUi(g,  h ) 


As  with  zAKl,  it  is  clear  that  AMi(g,  g)  =  0.  Calculation  reveals  that 
dgAMI(g,  h)\h=g  =  dhAm(g,  h)\h=g  =  0  and  that  the  correspond¬ 
ing  Riemannian  metric  on  M  is 
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